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Abstract 

We derive a (2+1) -dimensional multicomponent long wave-short wave resonance interaction 
(LSRI) system as the evolution equation for propagation of TV-dispersive waves in weak Kerr type 
nonlinear medium in the small amplitude limit, by following the lines of reference [14J for the two- 
component case. The mixed (bright-dark) type soliton solutions of a particular (2+l)-dimensional 
multicomponent LSRI system, deduced from the general multicomponent higher dimensional LSRI 
system, are obtained by applying the Hirota's bilinearization method. Particularly, we show that 
the solitons in the LSRI system with two short-wave components behave like scalar solitons. We 
demonstrate that the bright solitons in the TV-component LSRI system with N > 2 display inter- 
esting collision behaviour resulting in energy exchange among the bright solitons, whereas the dark 
solitons undergo standard elastic collision accompanied by a position shift and a phase-shift. Our 
analysis on the mixed bound solitons shows that the additional degree of freedom resulting due to 
the higher dimensional nature of the system results in a wide range of parameters for which the 
soliton collision can take place. 
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I. INTRODUCTION 



The nonlinear interaction of multiple waves results in several new physical processes 
[IH3]. It has been shown in the two layer fluid model that resonance between the long- 
wave component and short-wave component occurs when the phase velocity of the former 
matches the group velocity of the latter [4] . This is a ubiquitous phenomenon which appears 
in hydrodynamics [5] , bio-physics [6] , plasma physics [7] , and in nonlinear optical systems [8] . 
Though there are many studies on the long wave-short wave resonance interaction (LSRI) 
in one dimension [7HT2]. results are scarce for multicomponent higher dimensional LSRI 
system. In the context of nonlinear optics, the interaction of bright and small amplitude 
dark pulses in optical fiber is governed by the integrable Zakharov model [8] . 

The resonance interaction of long- wave with short-wave was first investigated by Benney 
for capillary-gravity waves in deep water [13]. In this case simple interaction equations 
cannot be obtained, because for deep water waves there is no wave in the long wavelength 
limit. However, simple interaction equations can be deduced in a stable stratified fluid for 
oblique propagation of long and short- waves [5]. The single component two-dimensional 
LSRI equation for a two-layer fluid model has been derived in Ref. [4] by using the multiple 
scale perturbation method and bright and dark type one- and two-soliton solutions have 
been reported. In Ref. [14J, Ohta et al derived the integrable two component analogue 
of the two-dimensional LSRI system as a governing equation for the interaction of the 
nonlinear dispersive waves by applying the reductive perturbation method. Very recently, 
the non-integrable three component Gross-Pitaevskii equations have been reduced to single 
component Yajima-Oikawa system by using multiple scale method [15]. In another recent 
work [16], the one-dimensional integrable two-component Zakharov- Yajima-Oikawa equation 
has been derived using multiple scale method and special bright-dark one-soliton solutions 
have been reported. 

In Ref. [17J, we have obtained general bright M-soliton solution, for arbitrary M, of 



the same integrable multicomponent (2 + 1)D-LSRI system (see Eq. (12)) considered in 
the present paper. Our earlier work [17J on bright soliton solutions of the multicomponent 
LSRI system shows that the role of long interfacial wave is to induce nonlinear interaction 
among the short-wave components resulting in non-trivial (shape-changing) collision be- 
haviour characterized by energy exchange among the two short-wave components. As far as 
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we know, for the first time in Ref. [TTJ we have identified the shape-changing/energy sharing 
collisions of bright solitons in a two dimensional integrable nonlinear system. This system 
will act as a potential candidate for realizing soliton collision-based computing and multi- 
state logic [T81420] . Now it is of interest to derive the general two dimensional TV-component 
equations describing the interaction of several short-wave packets with long- waves in a phys- 
ical set up and to look for other types of multicomponent soliton solutions. 

In recent years, much attention has been paid to investigate mixed (bright-dark) soli- 
ton dynamics in different dynamical systems including nonlinear optical systems and Bose- 
Einstein condensates [U |8j [T4HT61 12TH23] of coupled bright-dark solitons and to analyse their 
propagation properties and collision dynamics. In the present work, we derive the (2+1)- 
dimensional TV-component LSRI system governing the evolution of m short-waves and n 
long- waves with (m + n = N) in a nonlinear dispersive medium and reduce the system to an 
integrable system for a particular choice of the system parameters. Then by applying the 
elegant Hirota's direct method to the integrable multicomponent LSRI system which results 
for particular choice, we obtain the coupled bright-dark one- and two- soliton solutions. We 
will show that bright and dark parts of the mixed solitons in the two short-wave components 
case behave like scalar solitons whereas if we go for three and more short-wave components 
the multicomponent nature of the solitons will come into picture and one can observe inter- 
esting propagation and collision properties. It is straight-forward albeit lengthy procedure 
to extend the analysis to construct M-soliton solution, with arbitrary M. 

The present paper is organized as below. The general TV-component LSRI system is 
derived by applying the multiple scale perturbation method in the next section. In section 



III, bilinear equations for the integrable (2+l)-dimensional multicomponent LSRI system 



are given. Sections [TV] and [V] deal with the mixed one- and two-soliton solutions of the 
multicomponent LSRI system. The collision dynamics of the solitons are discussed in section 



VT| Discussion on the mixed soliton bound states is presented in section [VTTJ The final section 
is allotted for conclusion. 



II. THE MODEL 

To start with we obtain the general two-dimensional multicomponent evolution equation 
for the propagation of TV-dispersive waves in a Kerr type nonlinear medium (ex.: optical 
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fiber, photo-refractive medium) by generalizing the approach developed in refs. [HI [Hj for 
the two and three components case. The waves are assumed to obey the following weakly 
nonlinear dispersion relations 

ujj = LUj(Kj]L,j : |^| 2 , \A 2 \ 2 , \A N \ 2 ), j = 1, 2, 3, AT, (1) 

where Kj and Lj are the x and y components of the wave vector, Aj{= Aj{x,y,t)) and 
Uj are the complex amplitude and angular frequency of the j-th wave. The fundamental 
carrier wave is of the form e l ( K ox+h oy -ujot) ^ m0 st convenient way to derive the evolution 
equation for the amplitudes Aj's is to Taylor expand the angular frequencies cjj's around 
the x and y components of the wave vector of the carrier wave K and L , respectively, and 
the central frequency uj at \Aj\ = 0, as below: 

(ujj - u) ) = (^kJoAKj + (uj^o&Lj + ^(cj j -k j k j )oAK j 2 + ^(^. L J AL 2 

N 

+( Wj , KjLj )o(AK j )(AL j ) + £ KlA m p)o l^-l 2 + i = 1, 2, iV, (2) 

m=l 

where AKj = Kj — K , ALj = Lj — L , j = 1, 2, AT. In this Taylor expansion and in 
the following, the subscript '0' in the ( ) represents that the quantity appearing inside the 
bracket is evaluated at L = K , Lj = L , ujj = (jJq and \AA = 0. In Eq. (2), uijK- = — 



by replacing (cjj — cj ), AKj and ALj by the operators — Zt|, zj^ and ij^ 5 respectively and 
transforming to the moving co-ordinates x' = x— cJit, y' = y—Uj^t^ t' = £, with the assumption 
that beyond a particular component (say q th ) all the derivatives ^Ig 2 ;^ , j = g+2, iV, 
are same and so also the derivatives , i- e -> (^+i,K g+1 ) (= = (^<z+2,K g+2 ) = 

= (un,k n ) q = (say) and (^+i, Lg+1 ) = (^+2,L q+2 ) = ••• = (un,l n ) q = ^2 (say), and 
omitting the primes for simplicity of notation we get 



N 



iAjj + iVj X Aj^ x + ivjyAj^y + C-p^ Aj^ xx + Aj^ yy + O^Aj^y + ^ ^ i?p^|A/| Aj — 0, (3a) 

A/" 

^p,t + C^Ap^x + C^A p ^y + CftAp^y + ^ ^ 5^ | ^ | 2 A p = 0, (3b) 

j = l,2,...,g, p = g+l, g + 2,...,jV. 

Here the independent variables appearing in the suffixes after the comma denote partial 
derivatives with respect to that variables and the group velocities of the j th component along 



the x and y directions are Vj X = ((wj^^o—wi) an d v jy = ((^',1^)0—^2), respectively. Various 
quantities in the above equations (3) are defined as = ^— Uj,k ^ Kj ^ ) = ^— ( ^bh ^ ? 



(^i,|Af| 2 )o' ^ = 1,2,3, TV. 

We consider the case where the first g-components are in the anomalous dispersion region 
and the remaining (TV — g)-components are in the normal dispersion regime. Following Ref. 
[8] , the solutions of ^ are sought in the form 

Aj = ^(x,2/,t)e^, j = l,2,...,g, (4a) 
A p = (uo + a p (x, y, t ))S^ M *(*M , p = q + 1, g + 2, N, (4b) 

A/" TV 

where Sj = B^u^ A p = B^u^ and a p 's are assumed to take only small values. 

l=q+l l=q+l 

Substituting equations ^ in (|3^) and neglecting the higher order terms involving (a pi $ p ) 
and also their derivatives result in the equation 

K^j,t + v jx^j,x + v jy i Pj,y) + Cl^1pj,xx + ^2^j,yy + ^3^j,xy 

+ (i2 B i j) \^n^ + ( E 5 p °'H2tioap) ) ^ = 0, j = l,2,...,g. (5) 

In a similar manner, by incorporating ([4]) in ([3)3) and collecting the real and imaginary 
parts, we arrive at a set of two coupled equations. The resulting coupled equations can be 
grouped together to the following equation by differentiating the imaginary part equation 
twice with respect to and making use of the real part equation: 

Q 



a p,tt + C4^p,xxxx + C^CsCLp^xyy + C^C§a p ^ xxxy 



2 

3 \xx 



0=1 



N 

+2C 4 u 2 B^a hxx = 0, p = q + 1, q + 2, TV. (6) 

Equations ([5]) and ([6]) are general equations describing the two dimensional propagation of q 
waves in anomalous dispersion region and (TV — g)-waves in the normal dispersion region. For 
TV = 3 case with q = 2 the system has two-SW components and one long-wave component 
and coincides with the corresponding equations presented in Ref. [T4J . One can notice from 
the general form of equations ^ and ^ that for the same N{= 3) but with different q 
value, (say q — 1), there is another possibility which will have one SW component and two 



long- wave components, and ultimately the dynamics will be different from the q = 2 case. 
Also, this systematic generalization to TV component case is necessary to identify the way 
by which the additional wave components (modes) in the normal dispersion regime alter the 
governing equation for the three components case given in Ref. fT4] . 

To deduce an integrable equation associated with the combined systems ^ and ^ we 
choose all the £>^'s, l,p = q + 1, N, in Eq. |6) to be equal to a constant value (say —71, 
71 > 0). Then we get 



a p,tt + C^PjXXXX + C^C^CLrp^xyy + C^C§dp^ XXX y 



2 

3 \xx 



N 

-2C 4 7i^o Yl ai > xx = °' P = q + + 2 ' N ' ( 7 ) 

l=q+l 

In the following, we investigate the cumulative effect of the small amplitudes a p 's on the 
SW components by considering the superposition of these amplitudes involving only the 
sum of all the amplitudes and neglect all other combinations as they will be small due to 
the smallness of a p 's. Particularly, we add all the a-equations and define YliLq+i a i = L> By 
doing so we get 

/ N q 

Ltt + C\L XXXX + C/±C§L XX yy + C^C§L xxxy + C 4^0 I Y] Bj P 

\p=q+l j=l 

-2C 4 [N-q]u 2 0ll L xx = 0. (8) 

The dispersion relation for the linear excitation corresponding to the long- wave components 
is found as 



a 2 = c 2 k 2 



c z c z c z 



(9) 



where c 2 = 2C^u^ (N — q)^\. Note that the velocity of the linear excitation depend upon the 
number of components TV and increases as we increase the number of components (modes) 
in the normal dispersion regime. Thus our systematic generalization to TV-component case 
shows that by altering the number of components in the normal dispersion region one can 
change the velocity of the pulse. 

Next we apply the multiple scale approximation method to derive the two-dimensional 
multicomponent LSRI system as in Ref. [14]. We re-scale the variables t, x, y, a p and ipj as 

t" = et, x = Ve(x + ct), y" = ey, a p = e a p , ^ = e 3/4 S (j \ (10) 



where c is as defined above after the dispersion relation ([9]). Then the following set of 
equations results from Eq. rt8|) at the order e 5 / 2 



N 



2cL xt + C 4 u [ Y,B?\S U) L)=0, p = q + l,q + 2,...,N. (11a) 

\p=q+l j = l 



At the order of e 5 / 4 , we notice from Eq. (J5J) that all the group velocities of the short-wave 
components along the x direction are the same and their magnitudes are equal to the phase 
velocity of the long- wave component 'c' (i.e. Vj X — — c, j — 1, 2, g). This is the condition 
for resonant interaction between long- waves and short- waves. Equation ^ reduces to the 
following set of coupled equations at the order of e 7 / 4 after replacing vj x by c and rescaling 



of the variables as defined in equation(lO) 



(Sl j) + v 3y Sf) + d?S% + ( 2u jh BjpaA =0, j = 1, 2, q. (lib) 



In Eq. (11), after applying the transformations (10) the double primes in the new variables 



£ t', £ x' and £ y' are dropped, for convenience. 



Equation (11) is the multicomponent LSRI system in (2 + l)-dimensions which is non- 
integrable in general. By suitably choosing the constants Bj 's, Bp 's, C4, (7p , j = 
lj2,...,g,p = g+ l,...,iV, and 71, along with the assumption that there is no group velocity 
delay between the short-wave components, we arrive at the following (q + l)-component 
(2+l)-dimensional LSRI system for the 2-dimensional propagation of dispersive waves in 
weak Kerr-like nonlinear media, 

l (s¥ ) + S^)-S(£ + LS^ = 0, j = l,2,..., g , (12a) 

L t = ^J2\S U) \l (12b) 



In Eq. (12), the subscripts denote the partial derivatives with respect to those independent 
variables. As mentioned in the introduction, we have obtained more general bright M-soliton 
solution, with arbitrary M, of Eq. (fl2|) [17J. In reference [TTj, we have expressed the bright 



M-soliton solution of (12) in Gram determinant form and explicitly proved that the general 
multisoliton solution indeed satisfies the bilinear equations. We have also pointed out in 
the same work that for the two-SW components case (q = 2), the bright soliton solutions 
reported by Ohta et al in Ref. [14] follow as special cases of our general multi-soliton 



solutions [17] . As the (q + l)-component LSRI system (12) admits M-soliton solution, for 
arbitrary M [TT] , the system can be integrable [2U [25] . The study on the other integrable 
aspects of Eq. (fl2j) is under progress and will be published elsewhere. 



III. HIROTA'S BILINEARIZATION METHOD FOR THE (2 + l)D MULTICOM- 
PONENT LSRI SYSTEM 

Hirota's direct method is an elegant and powerful tool to construct the soliton solutions 
of a given integrable nonlinear evolution equation [24] . We look for the coupled bright-dark 



(mixed) soliton solutions of Eq. (12) in which the mixed soliton has m bright parts and n 
dark parts, such that m + n = q. To construct the mixed type soliton solutions, we perform 
the bilinearizing transformations, = S^ 171 ^ = L = — 2^(ln/), j = 1, 2, . . . , m, 
and / = 1, 2, . . . , n, (m + n = q), where gr^'s and /i^'s are arbitrary complex functions of 
x, y and t while / is a real function. The resulting bilinear equations are 

D 1 (gM.f)=0, j = l,2,...,m, (13a) 
D 1 (h®-f)=0, Z = l,2,...,n, (13b) 

/ m n \ 

D 2 (/./) = -2 ^^V^ + E^^* l ( 13c ) 



1 = 1 



where D\ = i(D t + D y ) — D x and D 2 = (D t D x — 2A), D-s are the standard Hirota's 
bilinear operators [24J, * stands for complex conjugation and A is a constant yet to be 
determined. One can have bright solitons for the choice A = in D 2 [HI E] and the 



bright soliton collision dynamics of system (12) has been discussed in Ref. [17]. However 



for non- vanishing 'A' values, the system can admit coupled bright-dark and dark-dark type 
soliton solutions. In this paper, we focus only on mixed (bright-dark) solitons corresponding 
to mixed type boundary conditions, that is, L x ' y ^^°° g(jn+i) x,y,t-^±oo cons j- an ^ 

j = l,2,...,m, / = l,2,...,n. 

This procedure can be very well applied to construct the dark-dark soliton solutions also. 
Here for convenience we consider the first 'm' short-wave components to be comprised of 
bright parts of the mixed solitons and the remaining n{= q — m) components to exhibit 
dark parts of the mixed solitons. To construct the mixed soliton solutions we expand the 
variables g^'s, /i^'s and / as power series expansions in a standard way [23j [24J. After 
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solving the resultant set of equations recursively we can obtain the explicit forms of g^'s, 
/i^'s and / and hence the multisoliton solutions can be constructed. 



IV. MULTICOMPONENT MIXED TYPE ONE-SOLITON SOLUTION 



The mixed one-soliton solution of (12) with m-bright and n-dark parts can be obtained 
by restricting the power series expansions as = X9i \ = (1 + X 2 ^)? / = 1 + X 2 /2 



and by solving the resulting equations, after their substitution into the bilinear equations 



(13) at various powers of \ recursively. The mixed one-soliton solution can be expressed in 



the following standard form, 



R 



= Ajk 1R sech ( VlR + ) e l ™, j = 1, 2, .., m, 



g(m+l) _ e i(Cl+<i>l+K) 



cos(0/) tank ( rj 1R + ^ ) + i sin ((/>/) 



-2k 1R sech rj 1R + 



R 



(14a) 

, / = l,2,...,n, (14b) 
(14c) 



The various quantities appearing in the above equations are given below: 

-l 



1 

4 

tan 



S |Q?,| v G 

-i ( hi ~ mi 



|pi| 2 cos 2 ( 



u 1R k 



IR 



, Aj 



km 

„2 i 2 



, VlR = k lR 



x + [ 2k u - ^ j y + 



k 



IR 



a\ e 2 

Wi + W) 



k 



IR 



Vu = hix - (k 1R - k n + uu)y + iv u t, Q = (ra, - b t )t + b t y + mix. 



(14d) 

(14e) 
(14f) 



In equations (14) and in the following the suffixes R and / of a particular quantity denote 
the real and imaginary parts of that quantity, respectively. Also a^'s, j = 1,2, ...,m, pi, 
uji(= Ui R + ioJu) and k\ = ki R + iku are complex parameters, while mi and fy, / = 1, 2, n, 

n 

are real parameters. The above solution is non-singular for the condition |p^| 2 cos 2 (^) > 



i=i 



^i R ki R . The amplitude (peak value) of the j bright part of the mixed soliton is Ajki R and 
that of (m + /)-th dark part of the mixed soliton is p\. The speed of the soliton is uj\ R /ki R . It 
can be noticed that both parts of the soliton have the same central position Rj2k\ R . But their 
phases are different. In fact, the phase of the dark component has two contributions, one 
from the background carrier wave and the other from (\)\. The quantity indeed determines 
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the darkness of the dark soliton. It is interesting to notice that the bright and dark parts of 
the mixed soliton of the LSRI system with more than two short-wave components display 
several interesting features in contrast to the case of just two-SW components as will be 
shown. To elucidate the understanding of such behaviour we present the explicit forms 
of one- soliton solutions for the two- and three- SW components and analyse them in the 
following subsections. For brevity, in the following we refer to mixed M-soliton solution 
with m-bright parts and n-dark parts as (rab-nd) mixed M soliton solution. 



A. Two-SW components (m = 1, n = 1) case 



This case admits only a simple type of bright-dark pair in which the bright part of mixed 
soliton appears in the first component and the dark part of the mixed soliton in the remaining 
component or vice- versa. The one-soliton solution for this case can be expressed as 



pi|W(<fr) - u 1R k 1R ) sech ( VlR + - ] e^ I+e \ 



R 



Pi e 



cos(©!) tanh [r] 1R + ^] + i sin^) 



where ^ 



In 



L = -2k\ R sech 2 ( r] 1R + 



2 v /\p 1 \2COS 2 (4>i)-w 1R k 1R 



R 



tan 



(15a) 
(15b) 
(15c) 



■ 1 (^),fl = tan- (|) 



, Ci = ( 



777,7 



bi)t + biy + mi x, and t]ir and rju are given in equations (|T4p) and (|l4f). 

The amplitude of the bright part ^^\/|pi| 2 cos 2 (0i) — oJiRku^J is independent of the pa- 



rameter but it is influenced significantly by the background carrier wave (pi). Such 
a mixed soliton at t = — 3 and y — — 1 is depicted in Fig. 1 for the parametric choice 



-3 + i, cji = 1 + 0.7i, m 1 = 1.4, p=l-i,b 1 = -0.2, = l+i. 



B. Three-SW components (q = 3) case 



Next we consider Eq. (12) with q = 3. For this case the mixed soliton can be split up 
into bright and dark parts among the three-SW components in two different ways. One 
corresponds to the (2b-ld) mixed soliton case where bright parts are in the and 
components while the dark part of mixed soliton appears in the component. The 
other possibility is a (lb-2d) mixed soliton case in which the bright part of the mixed 
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FIG. 1: Mixed one-soliton in two-component LSRI system. (IS^ 1 ^ 2 - solid curve, \S^\ 2 - dashed 
curve) 

soliton appears in the component while the dark parts are split among the remaining 
components and S^. 




1. (2b- Id) mixed one-soliton solution 



The one-soliton solution for this case where the bright parts appear in the and 
components while the third component comprises the dark part of the mixed soliton, 



can be written from (14) as 

= AjkiR sech ( r] 1R + ^ ) e 



R 



lT)u 



pi e 



«(Ci+0i+?O 



j = 1? 2, 



R 



cos(0i) tanh ( rj 1R + — ) + i sin(0i) 



(16a) 
(16b) 



where 
(/>i = tan 



-1 f feu -mi 



_ R 



1, 2, £i (m\—bi)t+biy+mix, R = In 



and T/ii? and 771/ are as defined in eqn. (14). L takes the same form as 
in eqn. (15:) with the above redefinition of R. Here one can observe that the a^-parameters 



4(| Pi |2COS 2 ((/)i)-u; 1 ^ih) 



appear explicitly in the amplitude of the bright soliton. The (2b-ld) one-soliton solution is 
characterized by twelve real parameters ag>, a« afl «<?, k 1R ,k u , oo 1R ,u u , P i R , P u, m u 
and &i and is restricted by the condition |pi| 2 cos 2 (0i) > oJiRkm for non-singular solutions. 
Such type of (2b-ld) bright one-soliton is shown in Fig. [2](a) at t = — 3 and ?/ = — 1 for the 
parameters k\ = 2 — 2z, o;i = — 1 — z, mi = 3, pi = 1 — z, &i —0.2, 1.8 + z, and 



a 



(2) 



1. One can also tune the intensity of bright parts without altering the depth of the 
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dark part of the mixed soliton by suitably choosing the parameter as can be evidenced 
from Fig. [^b) which is drawn for same parameter as that of Fig. [^a) except for = 1 + 
The soliton appearing in the long- wave component looks similar in both the cases and so we 
do not present it here. 

2.0 
1.5 

\S (lX3 f 1.0 
0.5 
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FIG. 2: Altering the intensity of bright soliton without affecting the dark soliton of (2b-ld) mixed 
one-soliton in three-component LSRI system by tuning parameter. (|S'^ 1 ^| 2 - solid curve, |5^ 2 ^| 2 
- dashed curve, |S( 3 )| 2 - dot-dashed curve 



2. (lb-2d) mixed one-soliton solution 



This case corresponds to the appearance of the bright part of the mixed soliton in the 
component while its dark part appears in the and components. The corresponding 
mixed one-soliton solution is 

S {1) = VlPil 2 cos 2 ( ( /> 1 ) + |p 2 | 2 cos 2 (0 2 ) -u 1R k 1R sech^ m + |)e^ + ^, (17a) 



i((l+(f>l+7r) 



= Pl e 



L = -2k 1R sech rj 1R 



cos(<pi) tank ( r] lR + ^ ) + i s'm(<f>,) 



R 



1 = 1,2, 



(17b) 
(17c) 



where R = In ( ( ^+ fe *) 2 ) -In (|pi| 2 cos 2 (0i) + |p 2 | 2 cos 2 (0 2 ) -uj lR k lR ), Q = (mf - b t )t + b t y + 

mix, (j)i = tan -1 { ^ u k ~^ 1 ^ / = 1,2, #1 = tan -1 (j^j, and rj 1R and r] U are as in Eq. (14). 
This solution is characterized by five complex parameters fci, pi, p 2 and uj\ and four 
real parameters bi and I = 1,2, along with the condition |pi| 2 cos 2 (0i) + |p 2 | 2 cos 2 (0 2 ) > 
UiRki R . It can be observed from the above solution that in contrast to the (2b-ld) case, 
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here the amplitudes of the bright and dark parts cannot be controlled by the a parameters. 
For illustrative purpose, in Fig. [3] we have shown the (lb-2d) mixed one-soliton solution for 
the parameters k\ = 1 — i, ui = —1 — i, mi = 2, m 2 = —2, p\ = 1 — i, p% = 1 + i, b± = —0.2, 
b 2 = 0.2, and a? = 0.2 - O.Oli at t = -3 and y = -1. 

2.0 
1.5 

ie(l,2,3)|2 

v I 1.0 
0.5 
0.0 

-10 -5 5 10 -io -5 5 10 

X X 

FIG. 3: (lb-2d) mixed one-soliton in three-component LSRI system. 




MULTICOMPONENT MIXED TYPE TWO-SOLITON SOLUTIONS 



It is a straightforward but lengthy procedure to construct the two-soliton solution. We 



obtain the mixed two-soliton solution of system (12) by restricting the power series expansion 



for gW% feW's and / as = X 9?+X 3 9%\ h® = (1+ X 2 h { 2 l) + X ^ ] ), f = l +x 2 / 2 +X 4 /4, 
j = 1, 2, ... ,m, I — 1, 2, . . . , n, and following the standard procedure [23J. The explicit form 
of (rab-nd) mixed two-soliton solution is given below. 



{j) e m +a ( 2 j) e n2 + e 



j = l,2,...,m,(18a) 



L 



1 r 



Pi 



/ = 1,2, ...,n, 



-2^(ln( D) ). 



(18b) 
(18c) 
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where 



JJ — l _|_ e m+V*+Ri _|_ e m+^2+ 5 o _|_ e ?72+?7t+^ e r]2+V2+ R 2 _|_ e Vi+V*+V2+V2+ R 3 



o (0 



Ci 



kjx — (ik 2 j + ujj)y + Ujt, j = 1, 2, 



(ki - k 2 )(a[ j) K2i - o^'ku) 
(h + kl)(k 2 + kl) 

(ki - ik) 



(k* + ih 



- jiij , i j j 1,2, e 



(k 2 - ki)(a^ 'k 12 - otf'n 22 ) 
(k 2 + k* 2 )(k 1 + k* 2 ) 

(h - ih){k 2 - ibi) 
[k\ + ibi)(k* 2 + ibi) 

r 3 _ \h - k 2 \ 2 (n U K 22 - /ci2«2i) 



(i), 



,Ri - 



/in, e R2 = fi 22 , e 5 ° = /ii2, e 5 ° = fi 21 , e 



jri P \^{h-imi){kl + imi) 



(h + kDlh + k^ih + k^) 



(h + K) 



(m, 



bi)t + biy + mix, i,p = 1,2; j = 1, 2, m; / = 1,2, ...,n. 



(18d) 
(18e) 

(18f) 

(18g) 

(18h) 

,(18i) 
(18j) 



We discuss below the two- and three- SW components cases to bring out certain interesting 
features of the multicomponent LSRI system with q > 2. 



A. Two-SW components case (m = 1. n = 1) 

For this case, the mixed two-soliton solution with its bright part appearing in the 
component and the dark part appearing in the component can be obtained by putting 
m = 1 and n — 1 in Eq. (18). The resulting (lb-Id) mixed two-soliton solution is 
characterized by seven complex parameters (a^ , , fci, &2, and pi) and two real 

parameters (6i and mi). This mixed two-soliton solution is restricted by the conditions 
\pi\ 2 {ki + k*) > (cji + (jO*)(ki — imi)(k* + mi), i,p = 1, 2, as in the case of the one-soliton 
solution. 



B. Three-SW components case (m = 2,n = 1) 

i. (26 — Id) mixed two-soliton solution 

First we write down the mixed two-soliton solution with its bright parts in the first two 
SW components and while the dark part appear in the third component S^ 3 \ For 
brevity, we mention the straight-forward procedure to write down the soliton solution from 



Eq. (18) instead of presenting the explicit cumbersome expressions. The explicit forms of 
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bright part of the mixed solitons appearing in and can be obtained by putting 



j = 1 and j = 2, respectively, in Eq. (183,). The dark part of the mixed soliton solution 



in the component results from Eq. (18 3) for 1 = 1. In a similar manner, the various 



quantities appearing in the (2b- Id) mixed two-soliton solution can also be obtained by 
choosing m = 2 and n = 1 in the expressions for fi ip and n ip . The two-soliton solution in 
this case is characterized by twenty real parameters with the conditions |pi| 2 (/^ + k*) > 
(ui + cj*)(ki - im 1 )(k p + ztoi), i,p = 1, 2. 

2. (lb — 2d) mixed two-soliton solution 

Another possible split up for the three SW components case is to have the bright part 
of the mixed two-solitons in one component (say S^) and the other two components (S^ 
and S^) comprise of dark parts. The obtained (lb-2d) mixed two-soliton solution can 



be deduced from Eq. (18) by putting m = 1 and n = 2. This mixed type two-soliton 

solution is characterized by sixteen real parameters and is restricted by the conditions 
2 \p l \ 2 (k i + fc*) 

> — ^— > (uji + cj*). i, p = 1, 2, for obtaining non-singular solutions. 

We wish to remark that our above analysis can be extended in a straight-forward way 
to construct three- as well as multi- soliton solutions. We have indeed obtained the mixed 
three-soliton solution, but we desist from presenting the solution here due to its cumbersome 
expression. Also, from the three-soliton solution we identify that the soliton collision is pair- 
wise and there is no multi-particle effects. Hence a detailed analysis of two-soliton collision 
is necessary as the higher order soliton interactions can be analyzed in terms of two-soliton 
collision. 



VI. SOLITON INTERACTION 

The multicomponent mixed type two-soliton solution presented in the preceding section 
contains all the information regarding the dynamics of two solitons in the (2 + l)D multicom- 
ponent LSRI system. To elucidate the understanding of the collision of mixed solitons, we 
present the detailed asymptotic analysis of the two and three component cases separately in 
this section. In particular, we study the interaction of solitons in the x — y plane. A similar 
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approach can also be very well applied to study the collision dynamics in the x — t plane. 



We choose the soliton parameters as k\ R > 0, k 2R > 0, ku > k 2 i, 



> 



> 



U2R 



, without loss of generality. For this choice, we find that for a fixed "t" if the soliton, 
say si, is localized along the straight line rji R = ki R x + (2ki R ku — uj\R)y + (jO\ R t ~ 0, then 
77 2j r will tend to ±oc as (x, y) —> ±oc. Similarly, if the soliton, say s 2) is localized along the 
straight line v\ 2R = k 2 RX + (2k 2R k 2I - uJ 2R )y + u 2R t ~ 0, then rj 1R -> ±oc. 



A. Two-SW components case 



A careful analysis of the asymptotic forms of the mixed two-soliton solution for the two 
SW components case shows that the intensities of the bright and dark parts of the solitons 
before and after collision remain unaltered. There occurs only a position-shift in the bright 
and dark parts of the two colliding solitons. The asymptotic forms of the solitons in these 
regions are given below. 
(i) Before Collision (x,y —> —oo): 
Soliton si 



S{ 1] - ~ ^{"sech (rn R + 



(2)- 



pie 



cos(^ 1 ' ) )tanh ( rji R + — ) + i sin(0^) 



Ri 



R 

L ~ -2k 1R sech z [r] 1R + ^- 



Soliton s 2 



S. 



(i) 



s. 



(2) 



~ A\ sech yq 2 R + 
L ~ -2k 2R sech 2 ( r] 2R + 



R3 — Ri 



2 

cos(<^ 2) ) tanh ( rj 2R 

R3 — R\ 



R3 — R\ 



i sin(0^) 



(19a) 
(19b) 
(19c) 



(19d) 
(19e) 
(19f) 
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(ii) After Collision (x,y — >> -\-oo): 
Soliton s\ 



S[ 1)+ ~ 4+sech + ^y^J e^, (19g) 
<?f )+ ~ pi e i (C 1 +^ 1) + 2 ^ ) ) cos(rf } ) tanh (77^ + + % sin ^) > (i 9 h) 

L ~ -2fc lfl sech 2 L fl + . (19i) 



Soliton 5 2 



^2 



(1)+ 



Af+sech^+^e'*', (19j) 



cos(0^) tanh ( r]2R + — ) + i sin(0^ 2 ^ 



(19k) 



L ~ -2k 2R sech 2 (v2r + y) ■ ( 191 ) 



where rfj,a^\ j = 1,2, -R3, <^n, and S21 are defined in Eq. (18), 0^' 

tan- 1 (^),/ = l,2. 

The amplitudes of the bright parts of the solitons s\ and s 2 before interaction (A{~ , ^4i _ ) 

(1) _ ^1 

and their amplitudes after interaction (A{ + , ^4? + ) are given by A\~ = , = 

( 5 11 _^1±^3 ( 5 21 _^2±^3 (1) ^2 

6 2 A\ + ^ — an d I • Here and in the following, the superscript 



2 ' 1 2 1 2 

(subscript) of A's denotes the soliton (component) number and — (+) represents the soliton 

before (after) collision. Upon substitution of the corresponding expressions for i?'s and 5's 

we find that the intensities of the bright parts of the solitons si and s 2 are same before and 

after collision. Similarly, we find the amplitudes of dark part of mixed solitons appearing 

in the component before and after interaction are same and are equal to p\. The 

two colliding solitons S\ and s 2 appearing in the bright component (S^) also experience a 

position-shift of opposite sense whose magnitude is given by |$| = | ^3-^2-^1 Additionally, 

the dark soliton si (s 2 ) experiences a phase-shift 20^ — tt (—20^ + tt). 

Thus the mixed solitons in the SW components undergo elastic collision. The solitons in 

the long- wave component also undergo elastic collision with mere position-shift of magnitude 

|$|. The elastic collision of solitons in the two component LSRI system is shown in Fig. [3] 

for the parametric choice k\ — 1 — i, k 2 = 1.5 + 0.75i, uj\ — 1 — 2i, uj 2 —0.75 + mi 2, 

pi = 2, bi = 7, ai 0.3 + z, and oip 0.05 — i at t — — 1. 
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FIG. 4: Elastic collision of (lb-Id) mixed solitons in two-component LSRI system. 



B. Three-SW components case 

This case can admit two types of soliton collisions as mentioned in section [V] B. The 
collision scenario depends upon the splitting of mixed solitons into bright and dark parts 
among the components and displays interesting dynamical behaviour. To illustrate this, we 
discuss the soliton collision for this (2b-ld) case in detail. 



1. (2b- Id) soliton collisions 



The asymptotic forms of solitons {s\ and s 2 ) before and after collision can be deduced 
from the exact two-soliton solutions (18) by putting m = 2 and n — 1. 
(i) Before Collision (x,y — »> —oo): 
Soliton si 



S{ j) ~ ~ Aj-sech (vir + y) elVU > i = 1 ' 2 
Sf ] ~ ~ pie 



cos(4>^) tanh ^rjm + ) + i sin(4>\ 

Ri 



Ri 



L ~ -2k 1R sech ( r) 1R + — 



Soliton S2 



M 3) - 



A 2 - sech T] 2 R + 



i?3 — Ri 



Pie 



*(Ci+^ 2) +2^ 1} ) 



L ~ -2k 2 R sech tj 2 r 



cos(0^) tanh ( tj 2 r 
i?3 — Ri 



3 = 1,2, 

i?3 — i?i 



i sin(0^ 2 ' ) ) 



(20a) 
(20b) 
(20c) 



(20d) 
(20e) 
(20f) 
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(ii) After Collision (x,y — >> -\-oo): 
Soliton s\ 



A - + sech ( t)ir + Rs R2 



J = 1,2, 



<(3)+ 



pie 



i(Ci+^ 1) +2^ 2) ) 



cos(<f>^) tanh ( t}ir + ^ 3 ^ 2 



-2k 1R sech 2 ( r) 1R + ^ 3 - ^ 2 



(20g) 
(20h) 
(20i) 



Soliton s 2 



°2 



Afsech + ^ 



^2/ 



J = 1,2, 



pie 



,i(Ci+4 2) +T) 



(^2) i -^2 

cos(0i ; ) tanh ^ 2j r + — 



i sin(0^) 



-2/c 2 r sech 2 [r]2R + ^ 



(20j) 
(20k) 
(201) 



The various quantities appearing in the above equations can be obtained from Eq. (18) by 
putting m = 2 and n — 1. Here, <^ = tan~ 1 ( fe ^^ ), / = 1,2. We find that the amplitudes 
of the bright parts of the two solitons Si and s 2 before and after interaction are related 
through the relation, 



A\ + = Tj>A l ~, j = 1,2, 1 = 1,2, 



(21a) 



U) 



where A)~ = \e~ R ^ 2 , A 



1 ft — Rl +i?3 

±e 1 lj 2 and the transition amplitudes T/'s are defined 



as 



h-h\ (k{ + k 2 



1 2 — 



^2 ^1 

h - k 2 



k\ k 2 J V /ci — I - /c< 



fci + A: 2 



1/2 



Ax-1 



1/2 



\/i — A1A2 



\/i — A1A2 

*^)A 1 



J = 1,2, 



(21b) 



(21c) 



in which Ai = ^ and A 2 = Here the /% p 's can be obtained from Eq. (18) by putting 



m = 2 and n — 1. This shows that the intensities of the bright parts of the mixed solitons 
before and after collision are different in general. The transition amplitudes T^'s, l,j = 1, 2, 
become unimodular only for the choice = \c^\/\c^\- On the other hand, the 
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intensities of the dark parts of the two solitons si and s 2 appearing in the third component 
remain unaltered after collision. Both bright and dark parts of the mixed solitons s\ and 
s 2 experience a position-shift of same magnitude but of different signs. The position-shift 
experienced by soliton si (and s 2 ) is $1 = ^3-^2-^1 (and $ 2 = _ $1), where R 2 and i? 3 



(2) 

are defined in Eq. (18) and the dark solitons si and s 2 experience phase-shifts 2(j)\ — tt and 



+ 7r, respectively. This shows that the bright components exhibit energy exchanging 



collisions characterized by an intensity redistribution (energy sharing) among the bright 
parts of the mixed solitons appearing in first two components and an amplitude dependent 
position-shift, whereas the dark parts of the two solitons undergo mere elastic collision 
accompanied by the same position shift as that of bright parts. Such a collision scenario is 
depicted in Fig. [5] for the parametric choice k\ — 2 — k 2 — 1.5 — U\ — 1 + i, uj 2 — — 1 + i : 
m 1 = 0.7, pi = 4, 61 = 2, = 1.2 + 0.2i, o£ ] = -1 + 2i, aj 2) = 0.25 + 0.25z, and 

= — 1 + ^ at t = 0. In Fig. 5J the intensity of soliton s\ is suppressed (enhanced) in the 
5(1) (g( 2 )) component and the reverse occurs for soliton s 2 - But the solitons appearing in 

and L components undergo elastic collision only. 




FIG. 5: Energy exchanging collision of (2b- Id) mixed two-solitons in three-component LSRI system. 
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2. (lb — 2d) soliton collisions 



Next we consider the collision scenario in the three-component LSRI system where the 
two colliding mixed solitons (say si and s 2 ) are comprised of one bright part and two dark 
parts and appear respectively in the and (S^ 2 \ S^) components. A careful asymptotic 



analysis of Eq. (18) with m = 1 and n = 2 shows that the amplitudes of the bright parts of 



the mixed solitons before and after collision are given by 



(A\~, At, A\\ At) 



a 



(i) 



2 / 



(1) 



(22) 



Substitution of the expressions for various quantities from Eq. (18) with m — 1 in Eq. (22) 



shows that the intensities of the bright parts of the two mixed solitons are same before and 
after interaction, i.e., |^4j + | 2 = |^4j | 2 , j = 1,2. Similarly, the amplitudes of the dark parts 
of the two mixed solitons before and after collision in the (S^) component are same 
and are equal to p\ (p 2 ). This clearly indicates that the intensities of the dark parts of 
the colliding mixed solitons are unaltered during collision. Thus for the (lb-2d) case from 
the above expressions we observe that both the bright and dark parts of the mixed solitons 
undergo standard elastic collision of solitons accompanied by position-shifts of magnitude 
|(i?3 — i?2 — i?i)/2|, where R 2 and i? 3 are defined in equation(18). The phase-shifts of 



the dark solitons si and s 2 are 2(ffi — tt and — 2(j>^ +7r, respectively. Such an elastic collision 
behaviour is shown in Fig. [6] for k\ 1 — 2i, k 2 1.5 + z, uj\ —1 — i ) uj 2 — 2 + i, vn\ —2, 
m 2 = 0.5, pi = 4, p 2 = 4, &i = 1, b 2 = 2, = 0.5 - z, and = 1.4 + i at t = -3. We 
do not present the plot for the long-wave component in Fig. [6] as it exhibits the collision 
process same as that of the component except for different amplitudes. 




x S2 15 





FIG. 6: Elastic collision of (lb-2d) mixed two-solitons in three-component LSRI system. 



Our analysis reveals the interesting fact that the energy exchanging collision of mixed 
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solitons can be realized only in the bright parts of the mixed solitons in three-component 
LSRI system and is possible only if the bright parts of the mixed solitons appear at least in 
two components. 

The above analysis can be extended straightforwardly to arbitrary NSW components 
case where the bright parts of the mixed soliton appear in m-components and the remaining 
(N — m) components admit dark parts. It can be shown that the shape-changing (energy 
exchanging) collision is possible only if the bright parts of the mixed solitons appear at least 
in two short-wave components, i.e. m > 2. 



VII. SOLITON BOUND STATES 



Soliton bound states are another interesting class of multisoliton solutions. Soliton bound 
states can be viewed as composite solitons moving with a common speed, two-soliton bound 



states in the (2 + 1)D two component LSRI system (12) with q = 2 can be obtained from 
Eq. (18) for the choice ^ = ^ and ku = k 2 i, m = 1, n = 1, and the corresponding 



solution reads as 
5« = 



5 (2) 



1 / h+ s 2i - h+ s n 

— (e 2 cosh(f) 2 R + iS 2 i) + e 2 cosh(r)i fl + iS n ) ) e mi 



Pie 



D 1 



R3 



A 1} 



e 2 cosh Ni 



; Q3I 



+ e 2 cosh(7V 2 + iQn) 



+e" 



Di= e 2 cosh r] 1R + t} 2 r + 



cosh(AT 3 + iQ 12 ) 



R1+R2 , / i?i 
+ e 2 cosh rim — t)2r H 



(23a) 



(23b) 



5 0+ S . . . 

+ e 2 cos(d /), 



where rj 1R = k 1R x + (2k u - %*)y + 



V2R 



(!) X 



u r 2L , Q 



11 - 
VlR + V2R " 



'111 



H Q 



12 



^12/ *3 



(1) 



o (1) 



VlR ~ V2R ■ 



VlR 

o w -o (1) 

^11R K <22R 



VlR + 7fe fl 



Vu - V21 ■ 



-o (1) 

112/ ^211 
2 



?7j/ = kjjx — (kj R — k 2 -j + ujji)y + cjj/t, j = 1,2, and all the other quantities found in the 



above expressions can be deduced from the corresponding quantities appearing in Eq. (18) 



by putting m = 1 and n — 1. The suffixes R and I appearing in the various quantities in 



Eq. (23) denote real and imaginary parts, respectively. 



The two-soliton bound state is shown in Fig. [7] (top panel) for the choice k\ = 1 + i, 
k 2 = 2 + z, uj r = -2 - z, cj 2 = -4 + 7z, mi 0.7, t -1, c^ 1} = 1 + z, 0.5 + z, p x = 4 
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and bi = 2. The bound state solitons display beating effects due to the oscillatory terms in 



Eq. (|23j). The beating effects can be suppressed completely by tuning the -parameters. 



This is shown in the bottom panel of Fig. [7j where the parameters are chosen same as 
that of the plots in the top panel except for q[, which is now fixed as 0.02. It should 
be noticed that though the a-parameters do not have any observable effects on the single 
soliton propagation in the two-SW components case as discussed in section |IV A[ they can 
display significant effects while considering bound state soliton propagation. 




FIG. 7: Mixed two-soliton bound states in the two-component LSRI system. Top panel: with 
beating effects, Bottom panel: suppression of beating by tuning the -parameters. 



Now it is of interest to investigate the influence of the cjj, j = 1,2, parameters on the 
bound state soliton dynamics which are arising due to the higher dimensional nature of the 
system. In this connection, we make u)i = u) 2 = and we choose ku = k 2 i. This will result 
in a bound state in the x — y plane which is stationary in time. Such a two-soliton bound 
state with breathing oscillations in the x — y plane is shown in the top panel of Fig. [8] for 
uji = uj 2 = 0, ki = 2 — z, k 2 = 1.5 — i, tyi\ — 0.7, pi = 4, b\ = 2, 1.2 + 0.2i, and 

a 2 — 1 + 2i at t — 0. But when uj\ and uj 2 become different and non-zero, they make 
the two solitons to undergo collision in the x — y plane. Thus due to the presence of the 
cjj -parameters and the higher dimensionality of the system there occurs a transition from 
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bound states to interacting solitons and is shown in the bottom panel of Fig. |8| This also 
shows that the presence of ujj parameters results in a wide range of parameters for which 
the soliton collision can take place. 
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FIG. 8: Mixed two-soliton bound state for uj\ — uj 2 — (top panel) and their transition to colliding 
solitons for uj\ — 1 + z, uj 2 — — 1 + z (bottom panel) in two-component LSRI system. Other 
parameters are given in the text. 



The dramatic change in the nature of soliton propagation in the x — y plane due to the 



presence of cjj-parameters resulting from the higher dimensional nature of the system (12) 
can exhibit additional features if we consider three SW components case. Soliton bound 
states m x — y plane for uj\ = u) 2 = 0, k\ = 2 — z, k 2 = 1.5 — z, mi — 0.7, p\ — 4, b\ — 2, 
= 1.2 + 0.2z, c4 1} = -1 + 2z, a{ } = 0.25 + 0.25z, and a£ 2) = -1 + z at t = is shown 
in Fig. [9j For the same parameters, with non-zero ujj values, there occurs head-on collision 
of solitons and this is depicted in Fig. [5] for uj\ = 1 + i and uj 2 = — 1 + z. Interestingly, the 
presence of second SW component now induces the fascinating collision involving energy 
exchange among the solitons in SW components. Also, the a-parameters can be tuned 
appropriately to suppress the beating effects of the bound soliton states. 
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FIG. 9: (2b-ld) mixed two-soliton bound states in three-component LSRI system for uj\ — UJ2 — 
and their transition to colliding solitons for ujj ^ is given in Fig. [5j 

VIII. CONCLUSION 

In this paper, we have derived the multicomponent LSRI system for the propagation of 
weak nonlinear dispersive waves in (2 + 1)D. Then we consider the integrable multicompo- 
nent LSRI system in (2+l)-dimensions and obtain mixed (bright-dark) one- and two-soliton 
solutions using Hirota's method. Our study shows that the bright and dark parts of the 
mixed solitons in the 2-SW components behave like scalar bright and dark solitons. But 
in three-SW components case different types of splitting of the mixed soliton into bright 
and dark parts are possible which makes their dynamics interesting. Study on the colli- 
sion dynamics of the mixed solitons shows that their collision in two-SW components case 
and in the three-SW components case with two dark parts and one bright part are elastic. 
However in the three-SW components case where the two colliding mixed solitons are split 
up into two bright parts and one dark part, the bright parts undergo energy exchanging 
collision characterized by intensity redistribution (energy sharing) and amplitude dependent 
position-shift. This collision process is also influenced by the soliton parameters of the dark 



25 



part. The dark parts of the mixed solitons undergo only elastic collision. 

Finally, we have considered the soliton bound states. To elucidate the understanding, 
we explicitly presented the two-soliton bound state expression for m — 1 and n — 1. In- 
terestingly, we find that the a-parameters which do not show any significant effect on the 
one-soliton propagation display interesting effects on the bound states. Particularly, they 
can be profitably used in suppressing the beating effects. Another important observation 
follows from our above study is that the presence of cjj -parameters can alter significantly 
the dynamics of solitons. Specifically, in the absence of cjj-parameters, which can result 
due to the higher dimensional nature of the system, there occurs bound state solitons in 
the x — y plane. But when the cj/s are brought into picture the solitons exhibit collision 
behaviour. Physically, this means that due to the presence of cjj-parameters the attractive 
force between the bound solitons vanishes and the solitons pass through each other. It has 
also been shown that for the three-SW components and one long-wave component case, for 
non-zero values of cjj, j = 1,2, one can have energy exchanging collision for the bright parts 
of the mixed solitons in the x — y plane whereas in the absence of u/s there occurs only 
stationary bound soliton in the x — y plane. 
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